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The physics beyond the Standard Model with parameters of the compressed spectrum is well 
motivated both in a theory side and with phenomenological reasons, especially related to dark 
matter phenomenology. In this letter, we propose a method to tag soft final state particles from 
a decaying process of a new particle in this parameter space. By taking a supersymmetric gluino 
search as an example, we demonstrate how the Large Hadron Collider experimental collaborations 
can improve a sensitivity in these non-trivial search regions. 


Introduction .—After discovering the Higgs boson at 
the Large Hadron Collider (LHC), the most important 
question is whether the scale of the physics beyond the 
Standard Model (BSM) is within the coverage of the 
LHC. As current LHC searches have been pushing away 
a possible energy scale of the BSM, most part of param¬ 
eter space around 0( 1) TeV mass scale in various models 
of BSM have been ruled out. But so far existing LHC 
analyses lose sensitivities in the compressed spectrum re¬ 
gion where the difference of masses among new parti¬ 
cles are negligible compared to their mass scale. Major 
difficulties in those analyses are from tagging soft par¬ 
ticles over Standard Model Quantum Chromodynamics 
(QCD) backgrounds. The parameter space of compresses 
spectrum is common in various new physics scenarios. 
In a supersymmetric framework (SUSY), the compressed 
spectrum is naturally predicted by various SUSY break¬ 
ing models [ , ! ]. For the extra dimension models like 
the Universal extra dimensional one (UED), the degen¬ 
erate spectrum in Kaluza Klein (KK) modes is usually 
predicted and only radiative correction can give limited 
mass difference between different KK particles [3]. From 
the phenomenological point of view, BSM with the com¬ 
pressed spectrum is favored. Dark matter annihilate pro¬ 
cesses in this parameter space case can easily satisfy a 
relic density compatible to current observations through 
the co-annihilation process as in the bino-wino or bino- 
gluino co-annihilation region [4, 5]. Due to the impor¬ 
tance of compressed spectrum, there have been lots of 
studies [(:]. Unlike previously suggested analyses, we pro¬ 
pose a method to tag soft jets from decays of BSM par¬ 
ticles. 

A fat-jet for compressed spectrum .— The event topol¬ 
ogy of our consideration is a three-body decay channel 
where a BSM particle A decays into two quarks and an¬ 
other BSM particle B , 

A —7 q + q' + B , (1) 
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with a case of compressed mass spectrum, 


Am = rriA — ms -C tua- (2) 


When the mass splitting is negligible compared to the 
mass of A, most of the energy from this decay flows into 
B making each quark q, q' too soft to be tagged as an 
isolated jet. This type of a decay is very common in 
BSM, some examples in SUSY with a gluino decaying 
into a neutralino or a heavier chargino/neutralino de¬ 
caying into the lighter neutralino when a mass splitting 
is less than W/Z/H bosons. We should note that cur¬ 
rent collider searches do not provide satisfactory results 
in this case, especially in search channels of multi-jet [ 
and of mono-jet+missing transverse energy (J£t) chan¬ 
nel [9] that barely provides limits on neutralino mass 
around 600 GeV. Thus we introduce a method to tag 
jets from A decay with suppressing the Standard Model 
backgrounds. To identify quarks in the above process, 
we propose a “fat-jet (FJ)” to capture quarks (q, q') as a 
single clustered object. The cone size of a fat-jet will be 
determined purely by the kinematics of a decaying pro¬ 
cess at the leading order QCD. A distance A R between 
two quarks would be, 


A 


1 M qq , 
y/z{\ - z) Pt( FJ) ’ 


( 3 ) 


where Pt(fj) is the total transverse momentum of ( q , q') 
(a fat-jet), z is a pp fraction of A that q obtains. In a 
case of a three body decay, due to a symmetric feature 
of (q 17 q'), z ~ 1/2. M qq / denotes an invariant mass of 
q, q'. With limited data, M qq * will be localised around 
the peak P qq i of an invariant mass distribution [10], 
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To select signal events over backgrounds, we require a cut 
on a transverse momentum pp of a fat-jet. With pprp j) > 
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FIG. 1: We show distributions of a distance A R between two quarks from a gluino decay (a) and an invariant mass M qq i 
(b) using a parton level Monte Carlo (MC) simulation. Here we choose Am = 40,100 GeV, solid lines with a requirement of 
Pt(fj) > Am and dotted lines without Pt(fj) requirement. 


Am, we have A R < V2 from eq. (3). Thus we set a 
radius R of a fat-jet to 1.5 to capture two quarks from 
A decay. In FIG. (1), we show a distance between two 
quarks from A decay and their invariant mass using a MC 
simulation at a parton level. We generate SUSY gluino 
pair productions and their decays into two quarks and 
a neutralino. The neutralino mass is fixed to 600 GeV 
and gluino mass varies according to the mass difference 
Am. As we see, we can capture most of two quarks from 
a gluino decay with a fat-jet radius R = 1.5 after a cut 
on a fat-jet, Pt(fj) > Am. Our approach has a benefit 
of statistical gain in a requirement on px of an initial 
state radiation (ISR) over mono-jet +$t searches. We 
only require ISR pt as large enough for A to be boosted 
resulting in Pt(fj) > Am. 

With introducing a large size jet, we have issues from 
soft QCD corruptions. A normal QCD jet (NLO) is ap¬ 
proximately to 

\J< Mj > NLO ~0.2 PJ R, (5) 

with a numerical factor 0.2 including colour charges [ ]. 

For an example with a mass splitting of Am = 100 GeV, 
a normal QCD jet attains a mass of •*/< Mj > ~ 30 GeV 
with Pt(fj) > Am. But as we show in FIG. (2), due to 
underlying events, background jets become more mas¬ 
sive. To remove these additional contributions on a jet 
mass and to sort out background QCD fat-jet 1 , we utilize 
a grooming technique and a jet substructure variable. 

Signal fat-jet v.s. QCD fat-jets. —We note that a sig¬ 
nal fat-jet is a two-prong jet while QCD fat-jet gets most 
of its energy from a single prong and obtains a volume 
by soft QCD radiations. A lot of studies have been 
performed to distinguish two-prong jets from the nor¬ 
mal jets [ - ] with focusing on large pt regions where 

preferences come from the motivation of boosted ob¬ 
jects. Here we show how one can adapt those techniques 


We use a term QCD fat-jet for the background QCD fat-jet. 


to the compressed spectrum in the region of moderate 
Pt- To show the difference between signal fat-jets and 
QCD fat-jets, we simulate supersymmetry gluino pair 
productions as in the previous section and Z+jet produc¬ 
tions for backgrounds. Both signals and backgrounds are 
generated by Madgraph_aMC@NLO [ ] with PYTHIA 
6.4.28 [16]. We apply ATLAS AUET2B [ 7,18] tune with 
PDF set CT6L] ]. For a detector simulation, we use 
Delphes v3 [20] interfaced with Fast Jet v3.0.6 [21] for a 
jet clustering with ATLAS detector parameters. Fat-jets 
are clustered according to Cambridge-Aachen(C/A) al¬ 
gorithm [22] with R = 1.5. 

To decontaminate a fat-jet from soft QCD corruptions, 
we implement the Mass Drop Tagger (MDT) [ ] out of 
various grooming techniques] 2, 23, 2 ]. We choose a 
parameters for MDT same as in BDRS Higgs tagger [ 
since our spectrums Am are about the similar order of 
the Higgs’ mass. In FIG. (2) we show jet mass distribu¬ 
tions of signal fat-jets and QCD fat-jets before and after 
passing MDT. There is a huge overlapping in jet mass 
distributions between signal fat-jets and QCD fat-jets be¬ 
fore MDT. However, after MDT procedure, the signal jet 
clearly shows a peak around 60 GeV, and the distribution 
is similar to our parton level analysis in FIG. (1) while 
QCD fat-jet becomes lighter. In short, our result shows 
that MDT method effectively removes soft-QCD corrup¬ 
tions in “non-boosted” pt region. In addition to this dif¬ 
ferent shape in jet mass distributions between the signal 
and the background, QCD fat-jets easily fail the MDT 
procedure since QCD fat-jet is characterized primarily 
by a single dense core of energy deposits surrounded by 
soft radiations that is not compatible to symmetric con¬ 
ditions. Thus MDT itself can act as a good analysis cut. 

To increase a tagging efficiency, we use a dimensionless 
parameter p = mj/ (p^i? 2 ) to impose a cut on fat-jets. 
This p parameter characterizes a splitting of two quarks 
inside a fat-jet. For a signal fat-jet with a radius R, p 
parameter becomes 
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FIG. 2: We plot fat-jet mass distributions of signal fat-jet and background QCD fat-jet before MDT (a) and after MDT (b) at 
a reconstructed level. For a signal process we use Am = mg — m^o = 100 GeV with m^o = 600 GeV. 



BP 

cuts 

Signal 

z + j 

m^o = 600 

Xi 

(etag ) 

(cfake) 

BPl 

10 < rriFj < 40 

0.75 

0.55 

Am — 40 

Additional p > 0.1 cut 

0.40 

0.18 

BP2 

40 < mFj < 100 

0.62 

0.22 

Am = 100 

Additional p > 0.1 cut 

0.42 

0.12 


TABLE I: Tagging/faking efficiencies on a fat-jet. We have 
two benchmark points (BP1, BP2) depending on a mass split¬ 
ting. Mass unit is a GeV. 


FIG. 3: p distributions after a MDT procedure at a recon¬ 
struction level. The study point of a signal process is the 
same as FIG. (2). 


from eq.(3). R qq > is the distance of the two quarks. For 
signal fat-jets, we can easily estimate a theoretical value 
of p ~ 0.2 with our choice ofpr(Fj) and R in the previous 
section. For QCD fat-jets, we take a quark-initiated jet 
as an example since most of jets in major background of 
Z+jet are quark-initiated. The leading order distribution 
of p from a jet consisting of two quarks [21 ]: 


p da 
a dp 


a s Cp 


©(p-S/*)ln 1 +0(y*-p) In — 
P V* 



The QCD fat-jet has a soft large splitting or hard small 
splitting which provides a small p value. In FIG. (3), a 
distribution of signal fat-jets clearly has a peak around 
0.2 as our estimation. The population of p in QCD fat- 
jet locates less than 0.1. This encourages us to impose a 
cut p > 0.1. The cut on p is independent on pt(fj) and 
provides an orthogonal constraint to the cut on fat-jet 
mass which is proportional to Pt(fj)- In Tab. (I) we show 
a tagging efficiency (etag) and a fake efficiency (ef a k e ) of 
tagging fat-jet. Signal fat-jet indeed has a much higher 
tagging efficiency compared to QCD fat-jet. 

Improvement on sensitivities .—We show a benefit from 
our method with a fat-jet by applying it to existing AT¬ 
LAS gluino search of the LHC 8TeV Run 1 [ ]. We select 


a signal region (2jm) in ATLAS search which is the most 
sensitive region in our benchmark points. An analysis on 
(2jm) requires at least two hard jets in final states and 
cuts are shown in Tab. (II). 

For our signal, we generate gluino pair production up 
to one jet matching using MLM scheme [26] to be con¬ 
sist with ATLAS analysis [ ]. For backgrounds, we simu¬ 
late three major backgrounds: Z(vv)- l-jets, W{tv)+ jets 
and tt+jets. For the Z(nP)+ jets, W(£^)+jets events, 
we generate Z/W production adding from 1 to 3 jets. 
For the tf+jets production, we generate tt with up to 
2 jets. For the hard QCD backgrounds, because of the 
requirement of the large missing energy and m e g, those 
QCD backgrounds are reduced to be less than 0.1% of 
the total backgrounds [7], so we can safely disregard its 
contribution to the total backgrounds. Additional jet in 
background events is considered through MLM match¬ 
ing procedure [ ]. Before applying ATLAS analysis, we 
require a fat-jet analysis with clustering particles accord¬ 
ing to C/A of R=1.5 and select a fat-jet with following 
criteria: 

(1) Apply a cut Pt(fj) > Am, on fat-jets (after MDT). 

(2) Choose the fat-jet with the largest p and mark this 
fat-jet as a candidate for a signal fat-jet. 

To show the effectiveness of above criteria, we define fol- 
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Signal region (2jm) 

cuts 

pT 

> 160 GeV 

Pr(ji) 

> 130 GeV 

Pr(ji) 

> 60 GeV 

minA <Xpt(., 1 i 2 , ( 3 )) ,#t) 

> 0.4 

/ y/ Ht 

> 15 GeV 1,/2 

m e ff 

> 1200 GeV 


TABLE II: ATLAS cut flows for signal region (2jm). 


lowing two efficiencies: 


ei 


-^signal 

Atotal 


£2 = 


-^candidate 


N f 


signal 


( 8 ) 


where TV-total is the number of generated signal events and 
-^signal denotes a number of events that contains a fat-jet 
from a gluino decay in a reconstructed fat-jet list after 
MDT. The Af canc jidate is the number of events where a fat- 
jet from a gluino decay is taken as the candidate through 
above criteria. Thus ei gives the fraction of events con¬ 
taining the signal fat-jets and £2 shows the efficiency of 
finding signal fat-jet. We find nearly half of the signal 
events contains the signal fat-jet and our method can 
reconstruct 60% — 70% signal fat-jets as in Tab. (III). 



AM = 40 GeV 

AM = 100 GeV 

ei 

0.64 

0.56 

€2 

0.58 

0.74 


TABLE III: Selection efficiency for the signal events. 


Methods 

Z+js 

W+js 

tt+js 

Total 

BP1 

a 

BP2 

a 

ATLAS 

430 

216 

47 

693 

74 

1.34 

57 

1.03 

With 

Am 

40 

48 

31 

7 

86 

27 

2.44 



FJ 

100 

21 

18 

11 

50 



17 

2.15 


TABLE IV: Expected number of events and corresponding 
significance of the LHC 8TeV with a luminosity of 20.3fb _1 . 
We refer readers to Tab. (II) for cuts of (2jm). The signal 
cross section for Am = 40 GeV is 0.83pb and for Am = 100 
GeV is 0.43pb. 

After preselecting events with mass window cut and p 
cut on a fat-jet as in Tab. (I), we recluster particles with 
anti-/c* [ ] of R = 0.4 according to the ATLAS analy¬ 
sis and apply cuts in Tab. (II). We compare the perfor¬ 
mance of our procedure of a fat-jet analysis with ATLAS 
analysis in Tab. (IV) 2 . The systematic uncertainty here 
is assumed to be 7% from in ATLAS (2jm) search [ ], 


2 The backgrounds without fat jet tagging are scaled to the back¬ 
ground events reported by [ 


800 

700 

vy 600 

o> 

O, 

?£ 500 

400 

300 

300 400 500 600 700 800 

nig [GeV] 

FIG. 4: We show an improvement from utilizing fat-jet tag¬ 
ging on search sensitivity of 2u in multi-jet signal region. 
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We would like to stress a point that we do not have any 
b— tagged jet veto which will reduce tt significantly [ ]. 

To see overall enhancement of the search sensitivity, 
we perform a complete scanning over ( rrig , o) as in 
FIG. (4). For a fat-jet analysis, we impose p > 0.1 cut 
and fat-jet mass window cut of 

-Am < tofj < A m for Am > 20 GeV, (9) 

5 GeV < m FJ < 20 GeV for Am < 20 GeV. (10) 

With the limited luminosity of the LHC Run 1, we gain 
only small increasement in the very degenerated case of 
Am < 20 GeV. This limit is from the jet qualification 
of pt > 20 GeV in pile-up removal [ 28 ] . With upcoming 
LHC runs, we will access the very degenerated region 
with a moderate help from ISR jet to have more fat-jet 
above the jet qualification. For the 14TeV LHC with 
a luminosity 30fb _ , we confirm the enhancement from 
fat-jet preselection as a significance of 2.4 compared to 
the result of conventional ATLAS analysis a ~ 1.64 in 
the studypoint of Am = 20 GeV with m^o = 700 GeV. 
In analyses using large size jets, pile-up may give effects 
on the filtered jet by smearing jet mass. But we note that 
pile-up would not affect our results significantly based on 
7TeV ATLAS analysis where the rate on the change of 
large size jet mass after MDT, = 0.1 ± 0.2 GeV up 

to 0(10) pile-ups [ ]. 

Conclusion. —We have proposed a method of cluster¬ 
ing soft jets from decaying process with a fat-jet in the 
scenario of the compressed spectrum. To identify signals 
over backgrounds, we apply a jet grooming method and a 
jet substructure variable to fat-jets. Tagging particles in 
signals will enable us to probe the properties (mass, spin 
and coupling structure) with upcoming high luminosity 
LHC [30-33]. There have been studies to use a tagged jet 
through jet substructure methods to see the possibility 
whether tagged jet maintains the parton level informa¬ 
tion [ ]. More detailed study in this direction should be 
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performed once the LHC observes the signal. 
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